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Abstract 

We consider a polyhedron with zero classical resistance, i.e., a poly- 
hedron invisible to an observer viewing only the paths of geometrical 
optics rays. The corresponding problem of scattering of plane waves 
by the polyhedron is studied. The quasiclassical approximation is 
obtained and justified in the case of impedance boundary conditions 
with a non zero absorbing part. It is shown that the total momen- 
tum transmitted to the obstacle vanishes when the frequency k goes 
to infinity, and that the total cross section oscillates at high frequen- 
cies. When the impedance Ao is real (i. e., there is no absorption), 
it is shown that there exists a sequence of frequencies k n such that 
the averages in the impedance of the total cross section over shrinking 
intervals around Ao go to zero as k n — y oo. 
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1 Introduction 



An interesting geometrical object was studied in a recent publication by 
Aleksenko and Plakhov [lj. This object O has the following property. Geo- 
metrical optical rays, coming from a particular direction and reflected twice 
from the boundary of O by the law of geometrical optics, continue to prop- 
agate parallel to each other in the same way as if the obstacle was absent. 
The object has zero classical total scattering cross section and appears in- 
visible to an observer on the basis of the theory of geometrical optical rays. 
Note that a phase shift may influence the "invisibility" of the obstacle. One 
should also note that optical ray considerations provide an approximation 
to the expected properties of the corresponding optical problem, when the 
obstacle is smooth and convex. These conditions do not hold for the object 
under consideration. A rigorous treatment of the problem has to be based 
on an investigation of the solutions of the wave equation. 

This paper concerns the study of the associated scattering problem for the 
reduced wave (i. e., Helmholtz) equation. High frequency asymptotics are 
obtained for the scattering of plane wave by the Aleksenko-Plakhov obstacle 

0. The rigorous justification will be obtained in the case of impedance 
boundary conditions with a non-zero absorption. It will also be obtained in 
a weaker sense for boundary conditions without absorption. It is important 
to note that the delay time is the same for all the rays meeting the obstacle, 

1. e., the phase shift A is a constant for all rays meeting O. In the absence of 
absorption, this implies that the obstacle O is almost invisible at a sequence 
of high frequencies k = k n = fc0 +^ 7rn ; n 00; where k° is determined 
by the boundary condition, and the invisibility effect disappears for other 
frequencies. 

Note that the scattering theory prohibits the existence of absolutely in- 
visible bodies, since a nontrivial outgoing solution of the Helmholtz equation 
can not have zero scattering amplitude. Our results show that the scattering 
data for one incident direction and a sequence of frequencies can be as small 
as we please for an obstacle of an arbitrary size. Another by-product of our 
results concerns the relation between the total scattering cross section and 
the geometrical cross section. For the obstacle under consideration, the total 
scattering cross section approaches to four times the geometrical cross sec- 
tion for some increasing sequence of frequencies (when the incident and the 
reflected waves are in-phase), while for convex obstacles the corresponding 
ratio approaches to two as the frequency increases. 

The simplicity of the model under consideration allows one to consider 
it as a candidate in competition with other models on cloaking (eg. |12j). 
Its obvious disadvantage, namely the incident direction is fixed, may be non- 
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essential in some applications. On the other hand, the corresponding cloaking 
object is easy to make, in contrast to other cloaking materials which present 
formidable technological difficulties. Indeed currently proposed cloaking ma- 
terials which lead to zero total scattering cross section for all incident direc- 
tions need to have an infinite mass [12] or to be highly non- isotropic [21 IT3] . 
The latter leads to huge difficulties in a practical realization of such a non- 
isotropic coating which are not attainable at present. 




Figure 1: Aleksenko-Plakhov object O. Translation distance \BD\ = 1. 

The Aleksenko-Plakhov obstacle C I 3 with Lipschitz boundary dO 
considered in this paper is presented in Fig. [TJ It is formed by an orthogonal 
translation of the 2D object shown in Fig. |2] The prolongation of this trans- 
lation is 1. In Fig. ^A'C'B' and ACB are equilateral triangles, where C,C 
are the mid points of segments A'B' and AB, respectively. The lateral sides 
A' A and B'B are slightly displaced in order to avoid a positive measure of 
the obstacle's surface on the tangent rays propagating along the z— axis. We 
suppose that A'B' = 1 that means that geometrical cross section of the body 
i. e.,the area of the projection of the body on the (x, y)-plane, equals to 1/2. 

The geometry of this obstacle is chosen in such a way that the geomet- 
rical optical rays propagating from the ^-direction, after a double reflection 
from the boundary of the obstacle, continue to move parallel to each other 
with a constant shift A of the phase. This creates the effect of invisibility for 
an observer who relies only on geometrical optical rays coming from above 
(we discuss this in more detail in the next section). In fact, the geometry 
of triangles ACB' and ACB can be chosen more generally, and other ob- 
stacles also can be constructed which also create the same effect of classical 
invisibility pQ. 

The goal of this paper is to study this invisibility effect using a rigorous 
analysis based on the investigation of the corresponding wave problem. It 
is well known that the justification of the geometrical optics approximation 
is a difficult task. This was done earlier in two cases: for smooth strictly 
convex obstacles QUI CHI El E] and for inhomogeneous media in the whole 
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space [13]. For example in the first complicated approximation of 

the field near the tangent rays is needed, and delicate local energy estimates 
are used. The latter are known for strictly convex or star-shaped obstacles. 
Some conditional results for non-convex obstacle can be found in [13] where 
the scattering amplitude is studied along the non-caustic directions. 

Let k > 0, the scattered field u is a solution of the Helmholtz equation 
satisfying the Sommerfeld radiation condition: 

Att(r) + k 2 u{r) = 0, r eQ = R 3 \0, (1) 



/ 

J\r\=R 



Oil \T) 

— iku(r) 



d\r\ 



dS = o(l), R^oo. (2) 



We assume an impedance boundary condition holds on dO, i.e. 

(J- + kX^(u + e ik ^) = 0, r = {x 1 y 1 z)edO ) Po = (0,0,1), (3) 

where > is a constant and n denotes the outer unit normal for O. 
There exists a unique solution in Hl oc (p) which satisfies all these conditions 
(eg [9]). Every solution u(r) of CO)-® has the following behavior at infinity 

e ik\r\ / 1 \ 

u(r) = ——UooiO) + o — , r ->■ oo, 9 = r/\r\ E S 2 , (4) 
\r\ \l r l/ 

where the function u^iO) = u 00 (^, A;) is called the scattering amplitude and 
the quantity 

2 / L, fn\\2. 



v(k) = ||moo||l 2 (52) = J 2 \ u oo(6)\ dfi(0) (5) 

is called the total cross section. Here dfi is the surface element of the unit 
sphere. 

The following observable is called the transport cross section. 



a T (k)= / (l-6- Po ))\ Uoo (6)\ 2 dfi(8) 
Js 2 

and equals the total momentum transmitted to the obstacle per unit time in 
large volume normalization. 

The following are the main results obtained in this paper. They are 
stronger if > (i. e., there is absorption of the energy at the boundary), 
and hold on average if = 0. 
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Theorem 1. Let^sX > 0. Then 

1) the transport cross section vanishes as k goes to infinity: 

lim oy = 0, 

k— >oo 

2) the total cross section has the following asymptotic behavior for large 

k: 

a(k) = ^\A 2 e ikA -l\ 2 + o(l), k^oc, A = (6) 
where A is defined in Fig® 

3) |woo(#)| 2 ~ ° r (k)5(po), k — y oo in the sense of distributions, i.e., for any 
<f G C(S 2 ), we have 

[ ip(6)\u 00 (6)\ 2 dS e = l-\-l + A 2 e tkA \\(po) + o(l), k -> oo. (7) 
Js* 2 

Remark. If A is real, then \A\ = 1. Thus if the second statement holds 
for real A then 

a(k) = - |-1 + e i(^+2arg(A))|2 + q ^ fc _^ ^ 

and 

2 arg(A) 2tt 

lim a(k n ) = for k n = h — n, n £ Z, SA = 0. 

71— >00 A A 

The almost invisibility of O manifests itself by the fact that the total cross 
section cr(k) can be made as small as we please by choosing a complex 
impedance A close enough to an arbitrary real A = Ao and then choosing 
k = k n large enough. 

It also follows from (E]) with real A that o~(k) — > 2 when k = ~ 2a ^s( j4 ) _|_ 
(2n+i)7r ^ j e ; a{k) approaches four times the geometrical cross section. 
The latter is the area of the shadow of the obstacle if the reflected rays are 
not taken into account, and this area equals 1/2. For convex obstacles, the 
ratio of the total to the geometrical cross section approaches two as k — > oo. 

In some cases we will denote the total cross section (j^J) by a\{k) in order 
to stress its dependence on the value of the parameter A in the boundary 
condition Q. 

Theorem 2. Let Ao be real and 

-2a rg (A(A„)) 2n 
A A 
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Then there exists a sequence of positive numbers e n — > 0, n — > oo, such that 



Consequence. Under conditions of Theorem [21 there exists a sequence 
of real numbers A„ — >■ Ao such that lim^oo a\ n {k n ) = 0. 

The paper is organized as follows. The eikonal geometrical optics approx- 
imation is constructed in the next section. It provides a basis to understand 
the validity of Theorems [H EJ It is also needed in order to construct the 
Kirchhoff approximation u° which is used to prove these results. Section 3 
outlines the proof of Theorem [TJ There we state the asymptotic properties of 
the Kirchhoff approximation in a bounded region and discuss the estimates 
needed to justify the approximation. In order to obtain the properties of 
u° we study first Kirchhoff approximation for the problem of scattering by a 
single polygon (one face of the obstacle O). This is done in section 4 with the 
main technical parts of the proof moved into Appendix 1. The asymptotic 
behavior in a bounded region of the Kirchhoff approximation for the obstacle 
O is obtained in section 5. Sections 6 and 7 provide the asymptotic behavior 
of the total and transport cross sections, respectively, when 3 A > 0. The 
case SA = is studied in section 8. 

2 Eikonal approximation 

We begin by recalling that the construction of the eikonal approximation to 
the solution of the scattering problem uses a Lagrangian manifold A formed 
by trajectories (r,p) = (r(t),p(t)) in the phase space R 3 x R 3 which cor- 
respond to the geometrical optics rays r = r(t). The components of the 
vector (r,p) are the position r = r(t) of a point along the ray at time t and 
the momentum p = p(t) (where p is the unit vector along the ray). The 
trajectories start at time t = at points (x ,yo,—a) of the plane z = —a 
with momentum p = (0, 0, 1). It is supposed that the obstacle O is situated 
in the half space z > —a. Location and momentum (r,p) of every ray, ex- 
cept those which meet an edge of O, are determined uniquely for all t £ R. 
Thus, (x ,y ,t) G R 3 are global coordinates on A. Let P r : A — > R 3 be the 
projection of A to the coordinate space R 3 and let B = P r dA C R 3 . 

In the case of the Neumann boundary condition, the eikonal geometrical 
optics approximation has the form: 





1/2 





: ikS{x ,y ,t) 



r G R 3 \£. (9) 



(x ,y ,t):r(x ,y ,t)=r 



6 



Here the action S is just the length of the trajectory, i.e. S(x, y,t) = t (since 
the exterior of the obstacle is homogeneous), the Jacobian \D(r)/D(x , yo,t)\ 
is the density of geometrical optics rays in ray tubes, and the summation is 
taken over all the points of A having r as projection to M 3 . 

Function ^ eik satisfies © outside B and is discontinuous on £>. Neverthe- 
less, we can use it to predict quasi-classic effects. Note that reflection from 
the plane faces of O does not change the density of the rays and so 



eik i 



(x ,yo,t):r(x ,y ,t)=r 



tkt 



r e R 3 \B. 



(10) 



In Figure [2] (where A = 1 in the case of Neumann or A = —1 in the case 
of Dirichlet boundary conditions) we see that the obstacle does not have a 
shadow: for example, the shadow zone of the right part of the obstacle is 
covered by the rays reflected from A' A" and B"B. 

The following observation on the geometrical structure is important: each 
trajectory started on the line A'B' between points A 1 and G or between point 
H and B' and ended on the line AB have the same length. It means that 
all the rays below the line AB have momentum po = (0, 0, 1), the action for 
the rays coming through GH is z, and the action for the rays which have 
collisions with the obstacle is A + z, where A = |CM"| + — |AM.|. So 

the rays below the line AB differ from the sets of rays in the absence of the 
obstacle only by the phase shift A on the rays coming through A'H U GH' . 
Hence, if kA is 



i multilple of 2ir 



the eikonal approximation <[)utside of some 
neighborhood of fhe obiftaate^eoir.gides witrHhe inciojgritewave J The scattering 




l^2gifc(z+A)l gifcz I yj2gifc(z+A) I 

Figure 2: Plane section of O. Base side jA'-B'l = 1. Eikonal approximation 

^eik is giyen above A , B , and below AB . A = + |^//_g| -\A'A\. 
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amplitude is not denned for the eikonal approximation, however the transport 
cross section for the eikonal approximation is defined by the following limit: 
lim^oo jj r|=fl (l - 9 ■ P oW lk (r) - e ikr ™\ 2 dS. 

Since all rays have forward direction after collisions with the obstacle, the 
transport cross section for the eikonal approximation is zero. The arguments 
above concern the case of the Neumann boundary condition (A = 0). When 
the impedance A is arbitrary, the terms in the eikonal approximation OH]) after 
each collision are multiplied by the factor 

ir^+A 
in a — A 

where n a = n-a is the cosine of the angle between the normal to the surface at 
the point of the incidence and the direction a of the ray before the incidence. 
In the case of our obstacle O all incidences occur with the same n a = —1/2. 

Define A(X) := |r§p then the eikonal approximation in the case of the 
impedance boundary condition can be written as 

ty eik (r) = [A(A)]" (xo ' 2/0 '* ) e ifet , r G R 3 \B. (12) 

(x ,y ,t):r(x ,y ,t)=r 

where n(xo,yo,t) denotes the number of collisions of the ray with initial 
coordinates (xo, yo) with dO before time t. Thus the field behind the triangles 
(after two collisions) is y 4 2 e lfc ( z + A ) ( see Fig. [2J). 

Note that the factor ffTTT) is chosen by the requirement for v|/ eifc (r) to 
satisfy the impedance boundary condition, i.e. 

^- + k\ ] m eik (r) = 0, re dO. (13) 
on J 

3 Outline of the proof of Theorem H 

The main technical difficulties of the present paper concern the construction 
of an appropriate asymptotic approximation for the scattered field. We will 
not work directly with ^> eik which has jumps and does not satisfy equation 
(TO . Instead we use the Kirchhoff approximation u° to the solution u in order 
to prove Theorem [2j This necessitates the need to justify the validity of the 
Kirchhoff approximation and to study its asymptotic behavior. 

Recall that the Green formula allows one to represent field u in terms of 
a surface integral involving u\do, f^lso- If ^O were smooth, the Kirchhoff 
approximation would be given by the same Green formula with u and its 
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derivative in the integrands are replaced by ?/> e * fc (r) := ^> eik (r) — e tkr ' Po and 
its derivative, respectively, i.e., 

i r a p ik\r-q\ f) p ik\r-q\ 

4vr J ac , <9n q \r - q\ dn q \r - q\ 



The integral above with e lkr ' Po instead of if} is equal to zero. Thus 

r f) P ik\r-q\ f) p ik\r-q\ 

/ ho-C^te)), +nr*( q ) MS q . 

Jdo 9n q \r - q\ dn q \r - q\ 



u m 

47T 

The exact definition of u°, in our case where dO is non-smooth, differs only by 
the introduction of a cut-off function 77 = T](k, q) in the integrand above which 
makes the integrand smoother and provides important uniform estimates of 
u°. Thus 

1 r a p ik\r-q\ X p ik\r-q\ 

A-k J do on q \r — q\ on q \r — q\ 

Here r] is a C°°-function defined on the faces of the polyhedron dO which 
vanishes in a /c -<5 -neighborhood of the edges of dO, equals one outside of a 
2A; _5 -neighborhood of the edges, and is such that |D m ?7| < C m k 5m , k > 1, 
for each partial derivative D m of order m on faces of DO. One can choose 
5 arbitrarily in the interval < 5 < 1/2. This choice of 5 allows one to 
apply the stationary phase method when an amplitude contains the factor 
7]. However, some of the arguments at the end of section 7 become simpler 
when 5 is small enough. Indeed it is sufficient to fix 5 = 1/4 and so 

\D m 7]\ < C m k m/ \ k > 1. (15) 

The function u° satisfies (111), Furthermore, it will be proved in section 
5 that the following statement holds 



Lemma 1. 

/ Ft \ 

= o(k), k — >■ 00. (16) 



-^ + k\) (u°{r)-u{r)) 



r)0 



L 2 {dO) 



This Lemma together with an a priori estimate for the solutions of the 
problem ([I])-© with > obtained in |4], allows one to justify the approx- 
imations u ~ u°, Uoo ~ u^, k — > 00. All the theorems in [3] were formulated 
for quite smooth obstacles. However, the proof of the Theorem 2 of that 
paper is based only on an application of the Green formula, which is valid 
for piecewise smooth obstacles also. Due to the importance (and simplicity) 
of this theorem for our purposes, we prove a slightly improved version of this 
theorem here. 
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Theorem 3. Let the function v satisfy £7]], ^) and the boundary condition 
d 



Then 



+ k\)v = f, r = (x,y,z) G dO, 3A > 0. 

on 



1 1 dv 1 . . . . 

|Poo||L 2 (52) < ^-^=||/||L 2 (ao), M\L 2 (dO) + \\-j^g^\\L 2 (aO) < -j^\\j\\L 2 (dO)- 

(17) 

Proof. Taking the imaginary parts of both sides of the Green formula 

/ (Av + k 2 v)vdx = I v n vdS — ik / {v^dS — / (|Vw| 2 — k 2 \v \ 2 )dx, 
Jn JdVt J s 2 Jn 

we obtain S J m v n vdS — k J s2 |t> 00 | 2 <iS' = 0. Thus, for any c > 0, 

k%\ [ \v\ 2 dS+k [ Iv^dS < [ \fv\dS<ck%\[ \v\ 2 dS+— ]— [ |/| 
Jan Js 2 Jan Jan ^ck^sX J dn 

Choosing c = 1, or c = 1/2 we arrive to the first inequality (jl7|) or to the 
second inequality for v, respectively. The second inequality for v n follows 
from the boundary condition. The proof is complete. 
Lemma [1] and Theorem [3] imply 

||wSo _ u oo\\l 2 (s 2 ) = o(l), k -)> oo, > 0. (18) 
Evidently for any </? G C(S' 2 ) we have 

/ v {0)\u oo {6)\ 2 dS e = I <^)|^(#)| 2 ^+o(l), fc^oo, > 0. (19) 

JS 2 

Hence, to prove Theorem [1] we need only prove Lemma [T] and to analyze the 
far field behavior of u°. 

Lemma [T]and the far field behavior of u° will be derived from the following 
lemma on the near field behavior of the Kirchhoff approximation. To state 
the lemma we need to introduce the set £> which is formed by extended to 
infinity lateral boundaries of shadow and reflected regions. One can define 
B as follows. Referring to Fig. El let dO 1 be the upper part of dO which is 
struck by the incident wave, and let DO 2 be the lower part of dO located 
strictly below dO 1 (it is struck by the reflected wave). Their left and right 
halves (rectangles) will be denoted by dO l s ,i = 1,2; s G {l,r}. Every time a 
geometrical optics ray strikes (directly or after a reflection) an edge of one 
of the four faces dO\,i = 1,2; s G {l,r} of the obstacle O, we extend the 
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ray to infinity in both directions. That is in the direction of the incident ray 
and in the direction of rays reflected from the face. These extensions form 
B. Consider, for example the incident ray (GA") on Fig. [2j Then B contains 
half-infinite ray (A"B) starting at A" and half-infinite ray starting at B and 
propagating downwards. 

Lemma 2. The following asymptotics hold uniformly on any compact of dO 
which does not contain points of edges and on any compact set in M 3 \i3 : 

u°(r) = ^ eik {r) + 0{1/ Vk), \Vu°(r) - V^ eik (r)\ = 0{Vk), k -)• oo. 

(20) 

For any R > there exist constants Ci = Ci(R),i = 1,2 such that 

\u°(r) - ip eik {r)\ < Ci, |Vu°(r) - Vip eik {r)\ < kC 2 (21) 
when \r\ < R, k — > oo. 

Remark. Note that ^ eik has jumps on B. In all estimates involving ^ eik 
on OO or B we assume that one-sided limits of the corresponding functions 
are considered. 

Lemmas |2] and [T] will be proved in section 5. After that Theorem [1] follows 
from the following two well known formulas: 



Q 



^-(r)+ik[ — -6)u u 



e~ ik{e - r) dS, 9 e S 2 , (22) 



M^<*> = (23) 

where Q is a closed surface such that the bounded part of the space with the 
boundary Q contains O. Note, one can take Q = dO. 



4 Scattering by a polygon McR 3 

This section is devoted to a study of a simplified version of the problem 
under consideration. It concerns the scattering by a single polygon which 
represents a typical face of the polyhedron O. The obtained results will be 
used later when the polyhedron O is considered. 

Let P C M 3 be a plane in M 3 , and let a G S 2 be a unit vector which is 
transversal to P. We choose a unit normal vector n G S 2 to P in such a way 
that n a = cos 9 = (n, a) < 0. 

Let M be a convex domain in P bounded by a polygon dM. By the 
shadow zone 5(a) = 5(a, M) C 1R 3 we mean the domain which is inaccessible 
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for the direct ray propagating along the vector a assuming that P\M is 
transparent and that M reflects the ray, i.e., 5(a) = {r = r + ta G M 3 : r G 
M, t > 0}. We define the reflected zone to be the area R(a) = R(a, M) C IR 3 
covered by the reflected rays, i.e., i?(a) = {r = r + ta* G M 3 : r G M, t > 
0}, where a* = a*(M) = a — 2(a • n)n G 5 2 is the reflection direction. In 
future a will be always equal to p (see Fig. [3J or Pq (see Fig. HJ). 

Denote by D{r) = D a : m( t ) the single layer potential with density ■qe lk ( a ' r ) ■. 

D(r)= / V e ik{a - q) j rdS(q), r G M 3 , (24) 

where 77 is the function introduced in ()14|) . For any r G M, define r e = r + rae, 
£ > and let 

d d 

—D(r) = \im—D(r e ). (25) 
On e^o dn 

Denote by N{r) = N a M(r) the double layer potential with density rje lk ^ a ' r ^ : 

JV(r) = / V e lk{a - q) — j -dS(q), q^M, (26) 

and define its value on M similar to (125 p . i.e. 

N(r) = limMrA 7— iV(r) = lim — MrA r G M. (27) 

Finally define A = A(a, M) = |^ (see (HU). 

Consider scattering of e* fcr Q by the obstacle consisting of the surface M 
only, and define the eikonal approximation ^^ k M as in f[T2"j) . Thus (see Fig. [3]) 

f _ e ifcr - a , r6S(a,I) 
^ k M (r) := *f M (r) - e ifcr ' Q = \ Ae lk ^^ +t °\ r G R(a, M) 

{ 0, R 3 \{i?(a) U5(a)} 

where the constant t = t (a, M) is defined by the relation 

e ik{(a* -r)+t ] _ £ ik(a-r) ^ r E M. (28) 

Since a* = a — 2(a ■ n)n, (J2"B"]1 implies that 

t = 2{a-n){n-r), r G M. (29) 
Note that n • (r± — r 2 ) = when r l5 r 2 G M. Thus, t is a constant. 
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The function ^ „ ^(r) has one-sided limiting values on both sides of M, 
the value of the function and of its normal derivative is zero on the shadow 
side of M. On the illuminated side M + of M: 

K%(r) = (A + l)e ikr - a , ^-^ k M ( r )=ik(-A + l)n a e lkr - a , r e M + . (30) 
Thus the Kirchhoff approximation ( I14j) in the case of M = O takes the form 

= $ a ,M = ^(A - l)n a D ajM + ~~~~~ N cl^m • (31) 

The following lemma plays a crucial role in the proof of the main result, 
theorem [TJ Informally it states that the Kirchhoff approximation in the case 
of O = M is close to ip^M := — e lkr ' a when k 1, and it justifies 

Fig. El 



<£ ~ 




Figure 3: The figure presents the main term of the asymptotics of the Kirch- 
hoff approximation <P = <P a ,M as k — > oo. Here R(a, M) is the reflected zone, 
S(a, M) is the shadow zone, and $ vanishes outside of R(a, M) |J S(a, M) 
when k — > oo. The main term of asymptotics coincides with 4>a M- ^ rigorous 
statement is given by Lemma [21 

Lemma 3. The following asymptotics as k — >■ oo hold uniformly on any 
compact set in M and on any compact set ofM. 3 \R(a, M) (J S(a, M) 

^,M(r) = <U0 + 0{1/Vk), |W„, w (r) - V<5*( r )| = O(Vfc). (32) 
For any i? > i/iere exzsi constants = Ci(R),i = 1,2 swc/i i/jai 

|^M(r) - < Ci, | W a , M (r) - V< fc M (r)| < kC 2 , (33) 

w/ien |r| < R, k — >• oo. 



13 



Remarks. 1) Note that nas j um P discontinuities on the boundaries 
of the reflected and shadow zones. In all estimates involving ^> eik on M or 
on the lateral sides of the boundaries of the reflected and shadow zones, we 
assume that one-sided limits of the corresponding functions are considered. 

2) Since ipa,M satisfies the impedance boundary condition (EJ), the esti- 
mates ( I3"2"j) . (13"3"j) imply that 

|-^,m + k\<P aM = -k(i(n ■ a) + \)e ik{a - r) + 0(Vk), k oo, 
on 

where the estimate of the remainder is uniform on any compact subset of M, 
and is bounded by C k on M. 

The proof of Lemma [3] is rather technical and is given in Appendix 1. 

5 Asymptotics of the Kirchhoff approxima- 
tion in a bounded region 

Proof of Lemma [2J Recall, see Fig. 121 that dOj, dO\ are the upper parts 
of the dO which are struck by the incident wave, and that dOf, dO^ are the 
lower parts of the dO which are struck by the reflected wave. 

Note that function ^> eik in (I14p vanishes on the part of the boundary 
of the obstacle which is not accessible for the rays, i.e., dO in (fl4l) can be 
replaced by the union of four polygons dOj, i = 1, 2, j = I, r. We split them 
in two pairs where the polygons in each pair are connected by rays. The first 
pair consists of Mi = 80] and Mi = dO^.. Then u° can be written in the 
form 

u° = £ + K, (34) 

where 

I r ft p ik\r-q\ 1 r X p ik\r-q\ 

£ = -t inr& ^ i i dS ^ / («) 1 \ dS » 

4tt J Ml u M2 dn q \r- q\ 4vr J Ml UMa dn q \r - q\ 

(35) 

and TZ is given by a similar expression with 80} |J 80^. in place of M\ [J M2 . 
We represent C as C = <&\ t i + <?2,r where the first term in the sum involves 
integration only over M\ and the second term involves integration only over 
Mi- Then $>\ \ coincides with the Kirchhoff approximation u° for a single 
polygon which is given by (1311) as studied in the previous section. In partic- 
ular, Fig. [3] and Lemma [3] with a = po are valid for $11. 

Function ^> eik in the region covered by the rays between M\ and M2 
has the form of a plane wave propagating in the direction Pq, i.e., $2,r is 
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<I> 



2.i 



I 

I 
I 
I 
I 



<P, r ~ 



Figure 4: The main term of asymptotics of <l>2,r 



also a Kirchhoff approximation for a single polygon. The only difference 
between <&\ r and $2,r is that in the second case the direction a of the incident 
wave is different (a = p$) and the incident wave contains a constant factor 
Ae ikt °. Thus, the main term of asymptotics of <£>2,r is given by Fig. HI This 
figure coincides with Fig. [3J with a = and extra factor Ae tkt ° added. An 
analogue of Lemma [3J is valid for <P 2 ,r stating that (£>2,r differs from the main 
term indicated in Fig. H]by 0(1/ \/k) on any compact outside of the lateral 
sides of reflected and shadow zones and by 0(1) on any compact set. The 
remainder is multiplied by k after differentiation. Thus, the main term of 
£ = + <3?2,r is indicated in Fig. [5j The shift of the phase t\ + t\ in Fig. [5] 
was discussed in section 2, where it was noted that this shift is equal to A 
(see Fig. [2]). Since 1Z has a similar asymptotic expansion and u° = C + 1Z, 
this proves the statements of Lemma [2J 

Proof of Lemma [TJ From ([2]) and (ITT1) it follows that 

+ kX ) ^ ~~ = °' re d °- ^ 

On the other hand, Lemma [2] shows that u° is close to ip eik (r) near DO. This 
leads to the statement of Lemma [TJ Indeed, let us fix an arbitrary e > and 
let dO e be such a small neighborhood of the edges on DO that the area \dO £ \ 

2 

of that neighborhood does not exceed 2(c 1 \x\+c 2 ) 2 ' wnere @h ^2 are constants 
defined in fj2Tj) . Then, for any k > 0, 

/ |(A + A;A)( M ° - ^' fc (r))| 2 ^ < (eA;) 2 /2. (37) 
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I 
I 

Figure 5: The main term of asymptotics of C 



Furthermore, (120]) implies that 



d 



+ kX)(u° - ^ eik {r))\ 2 dS <Ck< {ekf/2 (3* 

>dO\dO E vn 

if k is large enough. Statement of Lemma [T] follows immediately from fl36l 

(EHD. 



6 Asymptotic behavior of the total cross sec- 
tion 

Proof of the second statement of Theorem [TJ We begin with the 
evaluation of the total cross section of the Kirchhoff approximation u° using 
fT2"3"j) . Consider the boundary Q of a cube such that its faces are parallel 
to the coordinate planes and which contains the obstacle O. For any small 
5 > 0, we split Q in the following three parts: Q = Q1UQ2UQ3, where 
Qi is a 5-neighborhood of Qf]B (see definition of B in Lemma [2]), Q 2 is 
the orthogonal projection of the obstacle O into the low horizontal side of Q 
without points already included into Q±, and Q3 is the remaining part of Q. 

Note that area \Q 2 \ of Q 2 converges to 1/2 as 8 — > 0. Thus, from f[2"Tj) it 
follows that for any e > one can choose S such that 

1 

k 



. n du° . , „ £ , 1 . ^ , £ 

u°— < - and < - - Q 2 < -. 
or 4 2 4 



(39) 



With 5 fixed, one can choose ko such that 



Q 3 



<9r ~~ 4' 



(9?/° 



Q2 



<9r 



4 



(40) 
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where Z = \ — 1 + A 2 e lkA \ 2 . The latter relations follow from ( l2"Uj) . since (see 
Fig. [2D ip eik = ty eik - e ikz = on Q 3 , and 



^eik = A 2 e ik(z+A) _ e ikz on q 



1- 



Inequalities f l39]) . ( HUjl and formula ([23]) imply the validity of ([6]) for u°. 
Its validity for u follows from (fT8l) . The proof is complete. 



7 Asymptotics of the transport cross section 

The goal of this section is to show that the transport cross section for the 
obstacle under consideration vanishes as k — > oo if 3 A > 0. The last state- 
ment of theorem [T] follows immediately (see the last paragraph of this section) 
from the first two statements. As the first step, we need to obtain an integral 
representation for the scattering amplitude of the Kirchhoff approximation. 
Recall, that u° = C + 7Z (see flM])), where C is given by fl3"o"j) and 7Z has a 
similar form. Denote the scattering amplitude of the functions u°, C, 7Z by 
u^, Coo, TZoo, respectively (see @ where the scattering amplitude is defined 
for u). Recall that M 1 = 80} , M 2 = dO 2 , and M 2 is the shift of Mr, to be 
more precise, M 2 = M% + dp$, d > 0. 

Lemma 4. We have 

u° 0O = c oo + n oo , (4i) 

where 

= ^(9) [ ne^-oUds (42) 

47r J Mi 

with 

(0) = {A-l)n P0 -{A+l)n 6 + Aj kd{ - 1 -^ ((A - l)n po + (A + l)n e ) . (43) 

Function TZ^ is the reflection of with respect to the first argument, i. e 
Kootf) = £<*($), where 6 = 6(6) = (-0i,0 2 ,0 3 )- 

Proof. The validity of ( 14T]) follows from ( |34|) . so one needs only to find 
the terms in the right hand side of (l4Tj) . Scattering amplitude can be 
found from (I3"5"j) : 

£00 = ~ [ V^ eik (q))e- lke q dS q - ^ / r)^ ik (q)e- ik ^dS q 

47T JM1UM2 47T Jm!UM2 

= -— [ r](^-(^ etk (q))+tkne^ eik (q))e- lk9 - g dS q , 6 G S 2 . (44) 
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We need to insert here the values of the eikonal approximation \l/ e * fc (g) 
on Mi and M 2 . Its value and the value of the normal derivative on Mi are 
given by (130]) with a = po- In order to find the corresponding values on 
M 2 one needs to take into account the following fact which was discussed in 
section 5. The wave ty eik (r) shortly before the incident to M 2 has the form 
Ae lk ( p v' T+t °\ i.e., it differs from the original wave e lkpo ' r coming to Mi by the 
choice of the incident direction (pi instead of po) and an extra factor Ae lht °. 
Thus, \|> eifc (r) and its derivative on M2 are given by (130]) with a replaced by 
Po and the extra factor Ae tkt ° added in the right-hand sides of fl30]) . Hence, 
(jUJ) can be rewritten in the form 



In order to prove fj42|) . it remains only to rewrite the last integral above as 
an integral over Mi by using the substitution q — > q + dpi and the following 
two facts. Firstly, the normals n on Mi and M 2 have different direction, and 
therefore, n p *,rig on M 2 are equal to n Po , —rig on Mi, respectively. Secondly, 
since pi = po — 2(p -n)n and to is given by (|29|) where M = Mi, the exponent 
in the second integral above (with q replaced by q + dpi, q G Mi) can be 
rewritten as follows 



(PS ~0)-(q + dpi)) = (pS - 9) ■ q + d(l - p* ■ 9) 
= (po-9)-q-2(n-p )(n-q) + d(l-p* -9) = (p - 9) ■ q + d(l -p* -9) - t . 



The symmetry between and 7^00 can be observed from formula (1441) for 
Coo, and the corresponding formula for T^. 
The proof of Lemma H] is complete. 

Theorem 4. // > then the transport cross section of the Kirchhoff 
approximation u vanishes as k — >■ 00. To be precise, 



Proof. Denote the integral in (T42]) by D 00 (6 I ); it is equal to the scattering 
amplitude of the single layer f[24"j) . From Lemma H] it follows that it is enough 
to prove the statement of the theorem for the part of uq, i.e. Theorem H] 
will be proved as soon as we show that 






f3 po (9) = l-p -9, 3A>0. 

(45) 
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We do not need the exact form of the functions (3 po ,^, but only some 
estimates valid for those functions. Namely, 

0<(3 po (9)<C\d-Po\ 2 , OeS 2 , (46) 

|*(0)|<C, \y(6)\<C(\6-p*\ + k\6-p*\ 2 ), 9eS 2 , (47) 

In fact, is a real valued infinitely smooth function on S 2 with the 
minimum value at the point 9 — po- Thus, the validity of the quadratic 
estimate fT46|) is obvious. In order to justify fT4T|l . we represent \1/ in the form 
\]/(6i) = ty^O) + ^ 2 (9) where ^i(0) is given by ([43]) with the exponential 
factor in that formula omitted, and 

V 2 (0) = * (0) - # = A(e ikd ^- p o^ - 1) ((A - l)n w + (A + l)n e ) . 

Since n p * = —n po , function ^\(9) vanishes at 9 = p^. This function is 
infinitely smooth, complex- valued and /c-independent. Thus, |\l/i(0)| < C\9 — 
Pol- Further, 

|*2(0)| < C\e ik ^o {9) -1\< Ck\9-p*\ 2 . 



The latter inequality follows immediately from (I46|) . Hence, the validity of 
both flU and (jUD are justified. 

For any vector u E S 2 denote by u the orthogonal projection of u on the 
plane P containing M\. Vectors po^Po have the same projections p = p^. 
Let Bp~ (k~ J ) be the disk in P of the radius /c~ 7 , | < 7 < |, centered at po- 
Denote by fli, fl 2 small neighborhoods on S 2 of the points Po,Pq, respectively, 
whose projections on P coincide with B^ik^ 1 ). Then 

\9 - pol > *T 7 > when 9 eVl 3 = S 2 \(n 1 U Q 2 ). (48) 

Denote by A = A p the two-dimensional Laplace operator in the plane P. 
For any q, q G Mi and p = q — q <E P, we have 

^ e ifc(po-e)-g _ ^ e ifc[(po-0)-go+(p()-0)-p] _ _& 2 |po _ 0j 2 e *KPo-0)-?_ 

Hence, from (JUD and ffT5l) it follows that for any m 

1^00(0) = = — I nA m e ik(p °- e) - q dSa 

(ik) 2m \p -9\ 2m J Ml 

= ^— / (A m ri)e ik{po - e) - q dSa = 0(k~ 2m ^-^), 9 G 3 . 

(ifc) 2m |p - 0| 2m Af x 
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Since 7 < 3/4 and m > is arbitrary, it remains to show that 

k:=k 2 [ (3 P0 (9)\^(9)D oo (9)\ 2 dS = O(k- £ ), z = l,2, 3A > 0. (49) 

Let us prove (j49H for z = 2. The case i = 1 can be treated similarly (in 
fact, the latter case is simpler with the reference to (|46p instead of fT4Tj) ). 
We use Euclidean coordinates on P as local coordinates on f2 2 : 9 = 9(9). 

Since po,Po are transversal to P, the Jacobian J(0) = ^ = 1/ y^(l — I6 1 ! 2 ) is 
bounded when # G and 

|0-pSI < l^-Pol <c\e-pH 9 en,. 

From here, (T4T|) and the boundedness of /3(#) and |-Doo(#)| it follows that 
I 2 <Ck 2 [ \^(9)\ 2 dSg<C [ (k 2 \9 - pl\ 2 + k 4 \9 - pl\ 4 )dSe 

JB P0 (k-~<) JB P0 (k-~f) 

< 2nC / (jfeV + k 4 a 5 )da = 0(/T £ ), £ = 67 - 4 > 0. 
Jo 

The proof of Theorem H] is complete. 

Proof of the first and the last statements of Theorem [TJ The first 
statement of Theorem [T] follows immediately from f[T9"j) and Theorem |H Let 
us prove the last statement. It follows from the first statement that for any 
neighborhood Ss C S 2 of the point po, we have J S 2\ Sg \uoo\ 2 dS — > 0, /c — >• 00. 
Then from the second statement we obtain that 

/ \ Uoo \ 2 dS^ I |-1 + A 2 e ikA \ 2 , k->oo. 
Obviously, the latter two relations imply ([7]). 

8 Real impedance 

Proof of Theorem [21 Let us fix Ao G M. We will need the optical theorem 
for Aeffi: 

a x {k) = (4n/k)$ Uoo (Po), A G R. (50) 

Note that this theorem concerns the exact solution of the scattering prob- 
lem, i.e., it is valid for Uoo but not for the far field given by the Kirchhoff 
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approximation u^. Hence, asymptotic behavior of Qu^po) can not be de- 
rived from the known asymptotics of ||t4L>ll- We therefore calculate it inde- 
pendently. From lemma |4] we have 

«»(P0) = ^(Pd) / r](q)dS q =-n PQ (-l + A 2 e ikd ^^) [ V (q)dS q . 
^ Jm 71 Jm 

Note that d = \A'B"\ (see Fig. [2]) is the hypotenuse of the triangle A'B" H . 
Thus d(l - Po -p* ) = d- \B"H\ = A, and 

uKpo) = ^(-1 + A 2 e* kA )(l + o(l)), k oo. 

Hence, if k = k n defined by ([8]) and A is real, then on noting that A(X) does 
not depend on k and \A\ = 1 for real A, 

^(p ) = ^(-l + e 2 *^^^ k n ^oo, 3A = 0. 

(51) 

Our next goal is to estimate the difference f\(po) '■= M So(Po) — ^oo(Po) f° r 
complex A, 3 A > 0. First of all note that the relation ( 122]) for the scattering 
amplitude is valid also for u and v = u° — u. Thus, 

\f\(Po)\ < \\Vn\\La(dO) + 1 IH \l 2 (80)- 

Denote by A the set A = {A: |A — Aq| < 1, SA > 0}. Statements of Lemma 
[T] and Theorem [3] are valid uniformly in A G A. Together with the estimate 
above this implies that 

\fx( P o)\<^, k->oc, AG A, (52) 

where o(k) is uniform in A G A. 

Let e > 0. Note that min Ae(Ao _ eiAo+e )(A - (A - 2e))((A + 2e) - A) = 3e 2 . 
Thus, 

*m (A-(Ao-2 e ))((Ao + 2e)-A) 

<P(A) := > 1, A G [A — e, A + ej. (53) 

Now we are ready to estimate the average of the function (150]) . Using 
(153]) and the positiveness of the function (|50j) for real A and $(A) for A G 
(A - 2e, A + 2e), we get 



jr / 9fUoo(po)dA < — / ^(A)S Moo (p )rfA < — / <£(A)3 Moo (p )dA 

fc n JA -e K n J\ Q -. e fc ™ J\ -2e 
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A +2e i /-Ao+2e 

^(A)^(p )dA-- / <P(\)%f x ( Po )d\. (54) 

JA -2e ft n J A -2e 

From f[5~Tj) . the integrand of the first integral in the right hand side of 
fl54|) does not exceed Ce, i.e., the first term does not exceed Ce 2 , where the 
constant C is independent of e and k. 

Let us study the second part in (1541) . Since $(A) is a real- valued function 
for real A, we have 

-A +2e /-A +2£ 

$(\)%f x (6 )d\ = 3 / $(\)f x {0 o )d\. (55) 

'A -2e JXo-2e 

One can easily see that the Kirchhoff approximation u° defined by ffT4"|) 
is analytic in A. Then from (1221) it follows that u ^ is analytic in A. The 
function « M is also analytic in A, 3 A > 0, (see eg [3], [5]). Thus, f\(0 o ) is 
analytic when 3 A > 0, and the contour of integration in fl55|) can be replaced 
by the contour r = Ti U 1^ U defined in Fig. [6j 




A - 2e A A + 2e 

Figure 6: Contour T 



KA 



From here and 

"A +2e 



it follows that 



<P(\)fx(6 )d\ 



A -2e 



<p(\)fx(e )d\ 



< o(k) max^plrl 
~~ r 3A 1 1 



One can show that maxp 



WM < c 



Since |T| = 8e, we obtain the following 



relation which holds for the second term in the right-hand side of ( 15"4"|) . There 
exists e > such that 

"A +2e 

4>(\)f x (e )d\ 



i 



5 n := max — 

0<e<e k n 



A -2e 



— > 0, fc n — >■ OO. 



Now from (154")) and (|50|) we get that, for e e (0, Eq), 

1 <5 n 
77- / °x(k n )d\ <Ce + —, k n -+ 00. 

le JX -e £ 

The choice e n = ^f5^ l implies the statement of Theorem [2J The Conse- 
quence from the theorem follows if Ao is chosen to be the minimum value of 
o'x(kn) on the interval Aq — e n < A < Aq + e n . 
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9 Appendix I 



Here we prove Lemma [3j The proof is a direct consequence of formula fj3Tj) 
and lemmas [5] and [6] below concerning the asymptotic behavior of the simple 
and double layers as k — > oo, r is bounded. These asymptotics are different 
in the shadow zone, the reflected zone and in the complementary region. 
We show that on M and in shadow zone we have 

2-7T 

D(r) = e ik{a > r) + 0(£T 3 / 2 ), r e M U S(a), k ->■ oo. (56) 

If r belongs to the reflected zone, then 

D{r) = -^ e ^o e ik(a*-r) + 0^-3/2^ r G R ^ k ^ 00j (57) 

where constant to — to(a, M) defined from the continuity on M, i.e., 

e ikt 0e ik(a*-r) = e ik(*-r)^ r £ M. (58) 

If r is in the complement to those zones, then 

D(r) = 0(AT 3/2 ), k -)• oo, (59) 

However, these asymptotics are not uniform, and exact statement is as fol- 
lows. 

Lemma 5. 1. Expansion |56j) is valid uniformly on any compact subset of 
M and on any compact set strictly inside of the shadow zone. Expansion (57\ ) 
is valid uniformly on any compact set strictly inside of the reflected zone. 

2. Estimate / f55j) is valid uniformly on any compact set which does not 
have common points with the closure of the shadow or illuminated zones. 

3. The above expansions on compact sets strictly inside S(a) or R(a), or 
strictly outside of these zones can be differentiated any number of times with 
the remainders multiplied by k after each differentiation. 

4- Derivatives of D(r) of any order in directions tangential to M are 
continuous on M , and the expansion < f55|) on M can be differentiated any 
number of times in those directions (with the corresponding increase of the 
order of the remainders). The normal derivative of D(r) has a jump on M, 
and 

f) 

—D{r) = -2nr]e lkia > r) , r G M, (60) 
on 

where the left hand side is defined according to 
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5. The following uniform bound is valid on any bounded region of R 3 



\D(r)\<^^-, \r\<Ro,k^oo, (61) 
k 

\P m D(r)\ <C(R )k m -\ r $ M, \r\ < R , k ->■ oo, (62) 
where P m is an arbitrary partial derivative of order m. 

Proof. First, we discuss the smoothness of D(r) on M (obviously D G 
C°° when r ^ M). Since the integrand in (1241) has a weak singularity, D(r) 
is continuous on M. Further, let r £ M, q e M. For any unit vector v 
orthogonal to n (i.e. v is in the plane containing M), we have J^|r — q\ = 
— it-\t — q\ where -Jr, -rf- are derivatives in the direction of the vector v with 

OVq I ^ I OV ' OVq 

respect to r and q, respectively. Thus, after integration by parts we get 

— D(r)= / r]e ik{a ' q) ^ -dS q - / V e lk{a ' q) e iklr ~ ql — [-. AdS„ 

dv J M \r-q\ J M dv q \r - q\ 

V ,d 1 9 .y k i^)Mr-<i\) dS+H {v) =ik( a ,v)D(r) + H^\ (63) 



m \ r ~ q\ dv, 



q 

where = J M (^v) ^~^0~qT ^ dS q . Obviously is continuous, and so 
is Jj^Dfr). If the vector u is also orthogonal to n, then the same arguments 
imply 

8 2 

D(r) = -k 2 (a, u)(a, v)D(r) + ik(a, v)H u + ik{a, u)H v 



dudv 



Q2^ e ik((a,q)+\r-q\) 

du q dv q \r — q\ 



+ / a . L ., dS r (64) 



This provides continuity of second derivatives. The higher order derivatives 
can be treated similarly. 

The normal derivative of a single layer has a jump on the surface of 
integration. This fact and the limiting values of the normal derivatives can 
be found in the many textbooks on partial differential equations. In our case 
the limit (160]) of the normal derivative is proportional to the density and does 
not contain an integral term since the surface of integration is flat. Formula 
(|60|) can also be treated as a simple exercise in the distribution theory. 

Our next goal is to find the asymptotics of D(r), k — > oo, on a compact 
subset K C M. Note that the integrand in (I24p is not smooth in this case, 
and the stationary phase method cannot be applied. Let us introduce polar 
coordinates (p, (f>) on the plane P with the center at the point r and the 
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polar angle being zero along of the ray r whose direction coincides with the 
direction of the projection a of the vector a on P. If a is orthogonal to 
P (i.e. a — 0), then r can be chosen arbitrarily. Polygon dM in polar 
coordinates has the form p = p(0, r), < < 2tt, where p as a function of 
is continuous, 27r-periodic and analytic for all the values of except those 
which correspond to the vertices of dM. We write the integral D(r) in polar 
coordinates and integrate by parts in p. When r G K C M and k is large 
enough so that rj = 1 on K, we get 

/■2-7T pp(4>,r) 

D(r) = e ik{a ' r) \ \ r]e ikp{l+lslcos ^dpd(j), 
Jo Jo 

gife(a,r) i-2-k ;-27r rp(4>,r) f J,kp(<j>,r)(l+\a\ cos </>) 

?7p : , —dpdcf)), (65) 



zfc Jq l + |a|cos0 Jo Jo 1 + |a| cos 

Since vector a is transversal to P, \a\ < 1, and f n 27r -, . Jf — r = , 2n = ^ ZL , 

'II ' JO 1+|O!|COS0 ^/l— |S| 2 na 

formula (I65p implies 

1 27T /" 27r rP^' r ^ p ikp(ip,r)(l+\S\ cos </>) 

D(r) = -e ik ^ r \— - / / r? p - — dpd0), r <E K C M. 

ik n a Jo Jo 1 + M cos0 

(66) 

The interior integral, above after integration by parts m times, takes the 
form 



1 fPi'P^) P)m+l ikp((p,r)(l+\a\cos(j>) 

4^ m win o-i-i W i ^ = Q ( fc(1 " 3m)/4 )- 
Jo (ap) m+i (1 + |a| cos0) m+i 



The latter relation follows from (jToj) . Hence, for r £ K C M, fl66|) implies 

dsnD. 

The tangential derivatives of Z}(r) can be evaluated using (IBUj) . f l6"l|) and 
similar formulas for the derivatives of higher order. One only needs to note 
that functions similar to H^ v ' have order 0(k~°°), k — > oo, i.e., they decay 
at infinity faster than any power of k. This can be justified precisely as the 
similar statement above for the integral in ( 166]) . Thus, the asymptotics (|56|) 
on M can be differentiated in tangential directions. 

We now consider the asymptotics of D(r) and its derivatives on compact 
sets outside of the boundaries of the shadow and reflected zones. If r ^ M, 
the integrand in fj2^|) is smooth, and the stationary phase method can be 
applied. Let F = F(q) be the phase function in (pM|) . Since V g F = a — pEfr j 
the stationary phase points are defined by the equation 

r — q 

a — = cn, where c is a constant, q e M. (67) 

r — q\ 



25 



Both vectors on the left are unit vectors, and therefore, (IBTl) holds if and only 
if r — q = Ci« or r — q = c 2 a*, where Cj = |r — q\ > 0, q G M. This is possible 
only if r is in the shadow or reflected zone, respectively. The stationary point 
is on DM if and only if r belongs to the lateral boundary of these zones. If 
r belongs to a compact set whose intersection with the boundaries of those 
zones is empty, the stationary phase point is strictly inside M or outside of 
M. Thus, the stationary phase method implies the uniform validity of the 
expansions fl56l) . fl57|) . fl59|) for D and its derivatives on these compact sets. 

Since the amplitude factor rj(k,q)/\r — q\ in the integrand fl2^|) depends 
on k, we will justify the latter statement a little more rigorously. Consider a 
compact set J strictly inside of the reflected zone. Let J' be the projection of 
J into M parallel to the vector a*. Then J' is located strictly inside M, and 
the stationary phase point belongs to J' when r G M. Let q — q(r) G C°°(M) 
be an infinitely smooth function which is equal to one on J' and zero in some 
neighborhood of dM. Let k be so large that 77 = on the support of <j. We 
write D(r) as the sum D\{r) + D 2 (r) where Di,D 2 are given by fpMl) with 
additional factors 77, 1 — 77 in the integrand, respectively. The stationary phase 
method can be applied to Di, since the amplitude factor for that integral does 
not depend on k. In order to complete the proof of (I5T1) for r G J it remains 
to show that D 2 = 0(k~°°), k — >■ 00. If f(k,q) = rj(l — ?)/|r — g|, then 

B2 = £ /e ,^ = £ _^_ (VF . VF)ei ^ = £ v(^«. 

This integration by parts can be repeated as many times as we please. 
This combined with (FT3]) provides the estimate for D 2 . Thus, (15T|) for r G J 
is proved. Estimates for the compact sets in the shadow zone and outside of 
the shadow and reflected zones are treated similarly. 

It remains to prove the last statement of Lemma [5j We begin with the 
proof of flBTj) . The proof will be based on the same arguments which were 
used above in order to justify the asymptotics of D{r) on M. Let us assume 
first that r belongs to the reflected zone, i.e. r = q + ta*, q G M,t > 0. 
We introduce polar coordinates (p, 0) on the plane P with the origin at qo 
and the polar angle defined in the third paragraph of the proof (where polar 
coordinates were used). Then similar to the first equality ( 165|) . we have 

D{r) = e ik(a, qo ) / / I pdpd4, (68) 

Jo Jo g{p,t,(p) 

where / = p\a\ cos0+ ^Jt 2 — 2t|o ; |pcos0 + p 2 , g = ^/t 2 — 2t|o ; |pcos0 + p 2 . 

For any function u(p,t,4>), denote by u' its derivative with respect to 
p. Let h(p,t,4>) = f'(p,t,(j))/p. Since |S| < 1, one can easily check that 
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h{p,l,4>) e C°°, h(p,l,<j>) ^ and < c x p' 1 < h < c 2 p~\ p -)■ oo. 
Moreover, l/\hg\, \(hg)'\ < C < oo when t — 1, p > 0. Then from the ho- 
mogeneity of / and g with respect to (p, t) it follows that the same estimates 
(with the same constant C ) hold for all t > 0: 

t^-t, \(hg)'\ <C<oc, t,p>0. (69) 

Using integration by parts, we get: 

r2w np{<t>,qa) 

D{r) = e ma, qo ) / / / * , <f>)e ikf ^dpdcf> 

Jo Jo h{p,t,<f>)g{p,t,(f>) 

" ( } 6 A [ MP,t,0),( P ,t,0) e Jl - od0 

/•2tt fp(cf>,qo) „/ 1 

- (^)- 1 e 4fc(a ' 9o) / / (t L + ^(^ ) , )e ikfip ' t ^dpd(f). (70) 
Jo 7o % % 

We now split the last term into 1% + 12 with factors 7/ and 77 in the integrands, 

respectively. Formulas (1691) justify (1611) for the first term on the right hand 

side of (1701) and for J 2 . They also lead to the estimate I\ < Ck^ 1 f M \r]'\dS q 

which implies (16T|) for I\, since \r]'\ = 0{k 1 ^) and the support of r( has 

measure of order 0(A; _1 / 4 ). The proof of (16T1) when r is in the reflected zone 

is complete. 

The same arguments are valid if r in the same half space bounded by 
P, but does not belong to the closure of the reflected zone. Then q does 
not belong to M. The limits of integration in ( 168|) will be different, namely, 
(p E I, pi(4>,qo) < p < pi(0, go)- Here pi,p2 are values of p where the ray 
from the origin with the polar angle intersects DM, and / is the set of 
values of the angle for which such an intersection is not empty. This leads 
to (T701) where the first term on the right has to be omitted and the limits of 
integration in the second term have to be changed. Obviously, estimate (l6Tj) 
is still valid. In order to obtain (16T1) when r is in another half space, one 
needs only to replaced a* by a. 

Let us establish the validity of (16"2|) . We begin with the case of P m = D u , 
i.e. with the estimate of a tangential derivative of the first order. The 
desirable estimate follows from (|63|) : the first term in the right-hand side 
of (1631) satisfies (162]) due to (IBTl) . and the second term has an even more 
refined estimate. In fact, function H^ w ' has the same form as D(r), but with 
f] replaced by its derivative. Thus the arguments justifying (1621) for D(r) will 
provide the same estimate for H^ u ' with an extra factor k 1 / 4 in the right- 
hand side due to ()15|) . Tangential derivatives of D(r) of higher order can be 
obtained similarly. 
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We now estimate -j^D(r). Let r be the orthogonal projection of r onto 
the plane P, i.e. r = r + nt, r E P, G K. Then J^D(r) — I\ + ikl 2 , 
where 

I 1= [ r]—^—e ik{(a - q)+lr - ql) dS Q + [ ri-^—e^^-^dS,. 
Jm \r-q\ s J M \r-q\ 2 

The second factor differs from D(r) only by an extra factor in the 

integrand. This factor is bounded and homogeneous in (r, q) of zero order. 
Thus, the proof of flBTT) for D (see fl68l) - fT70l) ) can be repeated for I 2 , i.e. 
1-^2 1 < C/k, k — )• oo. Ii can be estimated very easily: if r% + t 2 < B%, (|r| is 
bounded), then by using polar coordinates centered at r$ we obtain 

i.e. (|62]) holds for £D(r). 

If P m contains differentiation of order m — 1 in tangential directions and 
the derivative of the first order in the normal to M direction, then (16"2"j) can be 
justified by a combination of the arguments used to prove (1621) for tangential 
derivatives and for ^D(r): one needs to start with fl63l) or similar formula 
for the tangential derivatives of higher order, then apply d/dn and repeat 
the arguments used to estimate -§^D. Finally, note that (A + k 2 )D(r) = 
when r ^ M. Thus 

— D(r) = (A M + k 2 )D(r), r <£ M, (71) 

where Am is the two-dimensional Laplacian on M. Hence P m D(r) for arbi- 
trary P m can be expressed through derivatives of D containing differentiation 
in the normal direction of at most first order. This proves fl62|) for arbitrary 
p 

The proof of Lemma [5] is complete. 

The double layer potential N(r) = N aiM (r) is given by fl2"E|) . f[2"Tj) . Ob- 
viously, since M is flat, N(r) = —-^-D{r). Thus, the next lemma about the 
properties of N(r) is a direct consequence of Lemma [5] and (T7T]) . 

Lemma 6. 1. If r 6 M and k — > oo, then 

N(r) = 2ue ik{a - r \ ^-N(r) = -2mk(n ■ a)e ifc(Q ' r) (l + OCl/Vk)), 
on 

where the estimate of the remainder is uniform on any compact subset of M 
and the remainder is bounded on M, i.e, \ 4^N{r)\ < Ck, r e M, k — > oo. 



28 



2. Let Ki, K-ii K 3 be compact sets strictly inside of the shadow zone or 
reflected zone, or strictly outside of these zones, respectively, and let t Q be 
defined by l[58\) . Then the following expansions are uniformly valid on these 
compact sets. 

N(r) = -27ce ik{a - r) + 0(1/ Vk), t G K\ C S(a), k ->• oo, (72) 
N(r) = 2ire iktQ e ikia *- r) + 0(1/ Vk), r e K 2 C R(a), k -> oo, (73) 
N(r) = 0(1/ Vk), reK 3 C R 3 \(R(a) U S(a)), k ->■ oo (74) 

These expansions on K{ can be differentiated any number of times with the 
remainders multiplied by k after each differentiation. 

3. A uniform bound of N(r) is valid in any bounded region: 

\N(r)\<C(R ), \VN(r)\<kC(Ro), \r\ < Rq, r $ M, k oo. (75) 
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